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We compute lattice dynamical properties of iron at the bcc-fcc phase transition using dynamical 
mean-field theory implemented with the frozen-phonon method. Electronic correlations are found 
to have a strong effect on the lattice stability of paramagnetic iron in the bcc phase. Our results for 
the structural phase stability and lattice dynamical properties of iron are in good agreement with 
experiment. 
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Elemental iron (Fe) is an exceptionally important 
metal for modern-day industry. For example, iron is 
the main ingredient of steel. Therefore a detailed under- 
standing of the properties of iron is essential for many 
technological applications. The remarkable macroscopic 
properties of iron are the result of a complex interplay 
between electronic and lattice degrees of freedom on the 
microscopic level. As a consequence iron exhibits a rich 
phase diagram with at least four allotropic forms. At 
ambient pressure the ground state is ferromagnetic and 
has a bcc crystal structure (a phase). Upon heating iron 
becomes paramagnetic at Tc ~ 1043 K. while retaining 
its bcc structure. Only when the temperature is further 
increased to T strU ct ~ 1185 K does the lattice transform 
to a fee structure (7 phase) [l[. Under pressure a iron 
makes a transition to a paramagnetic hep structure (e 
phase) at ~ 11 GPa. Even today, and in spite of long- 
term intensive research, the phase diagram of iron at high 
temperatures and pressures, as well as several key prop- 
erties, arc still poorly understood. 

State-of-the-art methods for the calculation of the elec- 
tronic structure provide a qualitatively correct descrip- 
tion of the equilibrium crystal structure and the lat- 
tice dynamical properties of the ferromagnetic a phase 
of iron. Various other properties of iron can also be 
understood on the basis of band structure calculations 
However, applications of these techniques to de- 
scribe, e.g., the a- j phase transition in iron, do not lead 
to satisfactory results. This is mainly due to the pres- 
ence of local moments above Tc which cannot be treated 
adequately by conventional band structure techniques. 
Namely, these methods predict a simultaneous transition 
of the structure and the magnetic state at the bcc-fcc 
phase transition while, in fact, the bcc-to-fcc phase tran- 
sition occurs only ~ 150 K above Tc. In addition, the 
elastic and dynamical stability of the bcc phase is found 
to depend sensitively on the value of the magnetization. 
For example, in the absence of a magnetization conven- 
tional methods find bcc iron to be mechanically unstable 
Clearly, an overall explanation of the properties of 
iron requires a formalism which can take into account the 



existence of local moments above Tc- 

The LDA/GGA+DMFT approach, a combination of 
the ab initio local density approximation (LDA) or gener- 
alized gradient approximation (GGA) of the density func- 
tional theory with dynamical mean- field theory (DMFT), 
allows one to determine the electronic and structural 
properties of materials with correlated electrons in both 
their paramagnetic and magnetically ordered states 0,0). 
Applications of LDA+DMFT have shown to provide a 
good quantitative description of the magnetization and 
the susceptibility of a iron as a function of the reduced 
temperature T/T c 0,0- 

It was found that the formation 
of local moments in the paramagnetic a phase is governed 
by the e g electrons which is accompanied by non-Fermi 
liquid behavior Q. This supports the results obtained 
with the s-d model for the a phase of iron Q. The 
LDA/GGA+DMFT approach implemented with plane- 
wave pseudopotentials [ijj II | has been recently em- 
ployed to compute the equilibrium crystal structure and 
phase stability of iron at the a-7 phase transition 12 1. 
The bcc-to-fcc phase transition was found to take place 
at ~ 1.3 Tc, i.e., well above the magnetic transition, 
in agreement with experiment. However, the lattice dy- 
namical properties of paramagnetic iron at the cv-7 phase 
transition, a problem posing a great theoretical and ex- 
perimental challenge, remained unexplored. 

In this Letter, we determine the structural phase sta- 
bility and lattice dynamics of paramagnetic iron at fi- 
nite temperatures by employing the LDA/GGA+DMFT 
computational scheme. The approach is implemented 
with plane-wave pseudopotentials which allows us to 
compute lattice transformation effects caused by elec- 
tronic correlations 1(3- 3 . We employ this ab initio 



computational scheme in combination with the method 
of frozen phonons [l3| to calculate the temperature de- 
pendent phonon dispersion relations and phonon spectra 
of paramagnetic iron at the bcc-fcc phase transition, a 
computation which was not possible up to now. 

We first compute the electronic structure of iron within 
the non-magnetic GGA b y e mploying the plane-wave 
pscudopotcntial approach [14[ and determine the equi- 
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librium lattice constant a for the bcc and fee structures. 
Overall, our results agree well with previous band struc- 
ture calculations @], e.g., we find a = 2.757 A for the bcc 
phase and a = 3.449 A for the fee phase of iron. It should 
be noted that these values are considerably smaller than 
those observed in the experiment Using the equi- 
librium lattice constants calculated above, we compute 
lattice dynamical properties of iron, by using the first- 
principles linear response method [1J] . In Fig. [1] (top) we 
show the results for the phonon dispersion curves and the 
corresponding phonon density of states of bcc Fe. The 
non-magnetic GGA finds the bcc lattice to be dynami- 
cally unstable with negative clastic constants Cn and C 
(see Table [TJ) , in accordance with previous calculations 
[3(. By contrast, the same method finds the fee lattice 
structure (with a = 3.449 A) to be mechanically stable 
(see Fig. [2]) and the calculated phonon frequencies to de- 
viate significantly from the experimental data. Indeed, 
the non-magnetic GGA finds a strong softening of the 
longitudinal [00£] mode at the X point by ~ 30%. Fur- 
thermore, all calculated elastic constants are 2-3 times 
larger than in the experiment (see Table [I]). Overall, 
calculations within the non-magnetic GGA find lattice 
dynamical properties which are in disagreement with ex- 
periment. Apparently, conventional band structure tech- 
niques cannot explain the experimentally observed struc- 
tural phase stability of paramagnetic iron at the bcc-fcc 
phase transition, since they do not describe electronic 
correlations adequately. 

To include the effect of electronic correlations, we com- 
pute an effective low-energy Hamiltonian for the par- 



tially filled Fe sd orbitals [llj, [18(. We employ results 



of band structure calculations for iron performed with 
the non-magnetic GGA to construct a basis of atomic- 
centered symmetry-constrained Wannier functions for 
the Fe sd orbitals [ll], [l8|- The corresponding first- 
principles multiband Hubbard Hamiltonian has the form 



H 



H, 



GGA 



- E 

2 V 

imm' ,oo' 



where hi 



; j TOCr Ci m(T , and c ima 



(cima) creates (de- 
stroys) an electron with spin a in a Wannier orbital m 
at site i. Here Hqga is the effective low-energy Hamilto- 
nian in the basis of Fe sd Wannier orbitals, and Hue is 
a double-counting correction which accounts for the elec- 
tronic correlations already described by the GGA. We 
use U = 1.8 eV and J = 0.9 eV in our calculations as ob- 
tained by previous theoretical and experimental estima- 
tions [19(. To solve the many-body Hamiltonian (1) we 
employ the DMFT together with quantum Monte Carlo 
(QMC) simulations with the Hirsch-Fye algorithm (20| . 

We now compute the lattice dynamics ofparamagnetic 
iron using the GGA+DMFT approach |4j, M in combina- 
tion with the method of frozen phonons [13| . The phonon 
frequencies are calculated by introducing a small set of 



displacements in the corresponding supercells of the equi- 
librium lattice which results in a total energy difference 
with respect to the undistorted structure. To this end we 
first perform a direct structural optimization and com- 
pute the equilibrium lattice constant of iron. We focus on 
the lattice dynamical properties of iron near the bec-to- 
fec phase transition. Namely, we perform our calculations 
at temperatures T = 1.2 Tc and 1.4 Tc, which are below 
and above the temperature T strU ct ~ 1-3 Tc where the 



structural phase transition occurs [12| . We calculate the 
total energy of the paramagnetic bcc and fee structures 
as functions of the volume and thereby determine the 
equilibrium lattice constants for the temperatures men- 
tioned above (see Table HJ. Our results for the equilib- 
rium lattice constants, which now include the effect of 
electronic correlations, agree well with experiment: for 
the bcc phase we find a = 2.883 A, which is only < 1 % 
smaller than the experimental value [l[, and for the fee 
phase a = 3.605 A, which is only < 2 % smaller than in 
experiment [l|. In the following we calculate the lattice 
dynamics of paramagnetic bcc and fee iron, respectively. 

Paramagnetic bcc iron: At T ~ 1.2 Tc we find the bcc 
phase to be energetically favorable, i.e., thermodynam- 
ically stable, with a difference in the total energy be- 
tween the bcc and fee phases of AE = Ef cc — E\, cc ~ 25 
meV/at. In Fig. Q] (bottom) we present our results for 
the phonon dispersion relations and phonon spectra. Wc 
note that these calculations are performed for the equi- 
librium volume (with a = 2.883 A) computed at this par- 
ticular temperature. To evaluate the phonon frequencies 
at arbitrary wave vector of the Brillouin zone we per- 
formed lattice dynamical calculations on the basis of a 
Born-von Karman model with interactions expanded up 
to the 5-th nearest-neighbor shell. The calculated phonon 
dispersions of the bcc phase of iron show the typical be- 
havior of a bcc metal with an effective Debyc temper- 
ature Od ~ 458 K. The phonon frequencies are overall 
positive which implies mechanical stability of the bcc lat- 
tice structure at T ~ 1.2 Tc, i.e., well above the Curie 
temperature. This result is quite different from that ob- 
tained with the non-magnetic GGA which finds the bcc 
lattice to be dynamically unstable (even for the equilib- 
rium lattice constant a = 2.883 A which almost coincides 
with experiment). Therefore this approximation cannot 
explain the bcc-to-fcc phase transition in paramagnetic 
iron. Most importantly, our results clearly demonstrate 
the crucial importance of electronic correlations to ex- 
plain both the thermodynamic and the lattice dynamical 
stability of the paramagnetic bcc phase of iron. 

Overall, the structural phase stability, equilibrium lat- 
tice constant, and phonon frequencies of bcc iron ob- 
tained by GGA+DMFT arc in remarkably good agree- 
ment with the experimental data which were taken at 
nearly the same reduced temperature T/Tc [3, EH. Nev- 
ertheless, we notice a weak anomaly in the transverse T\ 
acoustic mode along the [££0] direction, indicating that 
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FIG. 1: (color online) Phonon dispersion curves and corre- 
sponding phonon density of states of paramagnetic bcc Fe as 
calculated within the non-magnetic GGA (top) and DMFT 
(bottom). The results are compared with neutron inelastic 
scattering measurements at 1173 K [l5l ]. 



at T ~ 1.2 Tq the bcc phase may be close to an insta- 
bility. This result can be ascribed to a dynamical pre- 
cursor effect of the bcc-to-fcc phase transition. We note 
that a similar behavior of the T± [££0] phonon mode was 
found to occur in the S iron at high temperatures (22|. 
It appears that the temperature driven bcc-to-fcc phase 
transition in paramagnetic iron differs from the pressure 
driven bcc-to-hcp phase transition, where neutron stud- 
ies found no dynamical precursor effects 
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In addition, using our result for the phonon disper- 
sions, we compute elastic properties of paramagnetic 
iron. The elastic constants Cn, C12, and C44 (due to the 
cubic symmetry there are only three independent param- 
eters) are obtained from the estimates of the correspond- 
ing sound velocities along the [£00] and [££0] directions. 
Our results for the elastic constants arc summarized in 
Table U where they are compared with results from non- 
magnetic GGA and experimental data. The elastic con- 
stants obtained by GGA+DMFT are seen to agree well 
with the available experiments. 

Paramagnetic fee iron: Next wc calculate the lattice 
dynamical properties of the paramagnetic fee phase of 



TABLE I: Comparison between calculated and experimental 
elastic constants (in 10 12 dyn/cm 2 ) of iron. The calculated 
equilibrium and experimental lattice constants are presented 
in the last column. 



Method 


Phase 


T/Tc 


Cn 


C44 


C12 


C' 


a, A 


nm GGA 


bcc 




-0.15 


1.19 


4.92 


-1.87 


2.757 


DMFT 


bcc 


1.2 


2.30 


1.27 


1.57 


0.36 


2.883 


Exp.[l,i6] 


bcc 


1.1 


1.92 


1.24 


1.71 


0.10 


2.897 


nm GGA 


fee 




3.21 


1.97 


1.74 


0.73 


3.449 


DMFT 


fee 


1.4 


2.10 


1.38 


1.61 


0.25 


3.605 


Exp. [1, 12 


fee 


1.4 


1.54 


0.77 


1.22 


0.16 


3.662 



iron. Our calculations of the structural phase stabil- 
ity find that upon heating the fee phase is energetically 
favorable for T > T strU ct- The total energy difference 
between the bcc and fee phase is AE ~ —20 meV/at 
at T ~ 1.4 Tc- To prove the mechanical stability of 
the fee phase at this temperature we now compute the 
lattice dynamics of the fee phase of iron. Our results 
for the phonon dispersion relations and phonon spectra, 
which were obtained for the equilibrium lattice constant 
a = 3.605 A are shown in Fig. [5] (bottom). The ef- 
fective Debye temperature at T ~ 1.4 Tc is found to be 
Od ~ 349 K. The phonon frequencies are overall positive, 
implying mechanical stability of the fee lattice structure 
at T ~ 1.4 Tc- This qualitatively agrees with the results 
of non-magnetic GGA calculations which (for the GGA 
equilibrium volume) predict the fee lattice structure to 
be mechanically stable. However, the GGA energy for fee 
iron is higher than that for the close-packed hep struc- 
ture. By contrast, the GGA+DMFT calculations find the 
simultaneous thermodynamic and lattice dynamical sta- 
bility of the paramagnetic fee phase of iron at T ~ 1.4 Tc, 
in accordance with experiment. Our results for the struc- 
tural phase stability, equilibrium lattice constant, and 
phonon frequencies agree remarkably well with the avail- 
able experimental data taken at nearly the same reduced 
temperature T ~ 1.4 Tc ^ i s important to note 

that the application of the non-magnetic GGA to fee iron 
finds phonon frequencies which differ considerably from 
experiment. These findings clearly demonstrate the im- 
portance of electronic correlations for the lattice dynam- 
ical properties of fee iron. 

In conclusion, we employed the GGA+DMFT com- 
putational scheme to determine the equilibrium crystal 
structure and lattice dynamics of iron at the bcc-fcc 
phase transition. The calculated structural phase sta- 
bility and lattice dynamical properties of iron near the 
bcc-to-fcc phase transition are in overall good quantita- 
tive agreement with experiment. Most importantly, our 
calculations explain both the thermodynamic and the lat- 
tice dynamical stability of the paramagnetic bcc phase of 
iron below the bcc-fcc structural phase transition. In par- 
ticular, electronic correlations are found to be crucial for 
an explanation of the equilibrium crystal structure and 
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FIG. 2: (color online) Phonon dispersion curves and corre- 
sponding phonon density of states of paramagnetic fee Fe as 
calculated within the non-magnetic GGA (top) and DMFT 
(bottom). The results are compared with neutron inelastic 
scattering measurements at 1428 K [13]. 



lattice dynamics of iron. 
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